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Abstract 


For  a  system  of  an  M-level  atom  interacting  with  cavity  fields,  the 
statistical  properties  of  the  field  are  investigated  numerically.  The 
variation  of  photon  antibunching  and  probability  distribution  with  the  atomic 
level  number  and  initial  field  intensity  are  discussed  for  both  resonance  and 
off -resonance  cases. 
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I.  Int.v  .  Juccion 


There  are  two  aspects  of  the  problem  of  an  atom  interacting  with 

radiation  fields  in  a  cavity,  the  dynamical  behavior  of  the  atom  during  the 

interaction  process  and  the  evolution  of  the  coherence  and  statistical 

properties  of  the  radiation  field.  Although  a  large  amount  of  work  on  atomic 

properties  has  been  carried  out  .he  Jaynes -Cummings  (JC)  model, ^  relatively 

less  attention  has  been  paid  to  the  field  properties. 

Early  studies  of  the  first-order  correlation  function  in  the  jC  model 

revealed  that  the  disparity  between  <a^(t)a(t)>  and  <a^(t)><a(t)>  increases 

2 

with  time  if  the  field  is  initially  in  the  coherent  state.  Diagonal  elements 

3 

of  the  field  density  matrix  were  explicitly  calculated  later.  Only  recently 

were  effects  on  the  nonclassical  behavior  of  the  atom  and  the  field  due  to  the 

4 

cavity  loss  discussed.  The  upper  and  lower  bounds  of  the  Q-value  in  the  JC 
model  were  given  analytically,^  and  an  analysis  of  these  bounds  show  that  the 
larger  the  amplitude  of  the  initial  coherent  state,  the  less  sub-Poisson  (or 
super-Poisson)  the  photon  statistics  can  become. 

A  direct  generalization  of  the  JC  model  is  the  system  of  a  three- level 
atom  interacting  with  one-mode  or  two-mode  cavity  fields.  The  off -resonance 
behavior  of  the  cavity  fields  has  been  investigated  by  two  of  the  present 
authors  and  others. It  was  discovereo  '  r  Ref.  6  that  the  initially 
coherent  field  does  not  lose  its  coherence  right  after  the  interaction  with 
the  atom  takes  place  when  the  interaction  is  off  resonance.  Instead,  it  may 
recover  coherence  almost  periodically.  As  a  matter  of  fact,  the  field  can 
remain  coherent  for  a  long  time  if  the  detunings  are  sufficiently  far  away 
from  the  resonance.  The  photon  probability  distribution  and  fluctuation  are 
discussed  in  detail  in  Ref.  7,  where  it  is  found  that  for  a  one-mode  A-t3rpe 
atom,  sub-Poisson  and  super-Poisson  distributions  appear  alternatively  as  time 


develops,  and  for  a  two-mode  A-type  the  distribution  is  always  super-Poisson. 


Furthermore,  photons  involved  in  the  former  case  exhibit  antibunching 
phenomena  and  in  the  latter  case  exhibit  bunching  phenomena. 

The  photon  statistics  has  also  been  studied  in  another  case  in  which  an 
atom  with  one  common  upper  and  N-1  lower  levels  interacts  with  N-1  modes  of  a 

g 

cavity  field.  A  numerical  calculation  was  performed  for  resonance 
interaction  with  initially  one  mode  (called  the  pumping  mode)  of  coherent 
light  present.  Once  the  interaction  takes  place,  all  the  other  modes  (called 
signal  modes)  start  to  build  up.  Photons  in  these  signal  modes  are  found  to 
show  sub-Poisson  distribution  and  hence  exhibit  antibunching  all  the  time. 

Q 

As  the  observation  of  the  one-atom  maser  has  become  practical,  and 

since  the  atom  is  excited  to  high  levels  with  principal  quantum  number  in  the 

range  of  30  -  40,  it  is  inevitable  that  multiphoton  processes  are  involved. 

In  fact,  a  micromaser  based  on  degenerate  two-photon  transitions  is  already 

realized.  Theoretically,  however,  there  does  not  seem  to  have  been 

sufficient  attention  given  to  multiphoton  processes  in  the  literature  thus 

11  12 

far.  In  two  recent  papers,  ’  the  present  authors  considered  a  single  M- 
level  atom  interacting  with  cavity  fields  with  arbitrary  detunings.  Time 
evolution  of  atomic  level  occupation  probabilities,  mean  photon  number  and 
field  squeezing  have  been  investigated  in  detail.  Their  dependence  on  the 
atomic  level  number  as  well  as  on  the  initial  field  intensity  has  been 
discussed.  Here  in  this  paper,  we  discuss  the  photon  probability  distribution 
and  fluctuation  and  their  variations  with  various  parameters. 

An  outline  of  the  theory  is  given  in  Sec.  II.  Results  of  the  numerical 
calculation  are  presented  and  discussed  in  Sec.  Ill,  and  concluding  remarks 
are  given  in  Sec.  IV. 


II.  Theory 

The  general  formalism  of  the  theory  can  be  found  in  Ref.  11  (hereafter 
referred  to  as  I),  and  we  only  outline  the  essentials  here.  We  consider  a 
cascade  atom  with  M  energy  levels  as  shown  m  Fig.  1.  The  total  Hamiltonian 
is,  in  the  rotating  wave  approximation, 

M  M-1 

H  -  tifla^a  +  ^  ^  ^i^i+l^i^  +  h.c.  (1) 

i-1  i-1 

4-  * 

where  the  operator  a‘  creates  a  photon,  A^  creates  an  atom  in  the  i-th  level, 
and  A.  are  the  atom- field  coupling  constants.  Our  discussion  will  be  limited, 
for  simplicity,  to  processes  involving  only  a  one-photon  transition  between 
adjacent  atomic  levels. 

We  define  the  total  excitation  number  operator 

M 

N  -  a^a  +  ^ 

i-1 

which  is  the  sum  of  the  photon  number  and  atomic  excitation  number  operators. 

A  A 

It  can  readily  be  verified  that  [H,N]  -  0,  i.e.,  N  is  a  constant  of  motion. 

If  we  are  only  interested  in  the  particular  case  in  which  the  atom  is 
initially  in  its  highest  level,  then  we  need  only  consider  that  part  of  the 
Hilbert  space  involving  state  vectors  corresponding  to  N  >  M,  where  N  is  the 

A 

eigenvalue  of  the  operator  N.  This  is  a  consequence  of  the  fact  that  N  is  a 
conserved  quantity. 


Consider  the  state 


3 


|i,n>-(i>|n>  (3) 

in  which  there  are  n  photons  in  the  field  and  the  atom  is  in  the  i-th  state. 
The  subspace  corresponding  to  N  =■  M  +  n  is  spanned  by  the  state  vectors  lM,n>, 
|M-l,n+l>,  .  .  .  I  i  ,n+M-i>,  .  .  .  1 1  ,n+M- 1> ,  where  n  is  the  photon  number  in  the 
field  and  can  take  any  positive  integral  value.  Since  state  vectors  in 
different  subspaces  are  necessarily  orthogonal,  the  total  Hamiltonian  can  be 
diagonalized  in  every  subspace. 

An  arbitrary  state  in  the  subspace  corresponding  to  N  can  be  expressed 
as 


M 


1^  >  -  >  C.  , I i,n+M-i>  , 

l_  i.n+M-L* 

i-1 


(4) 


where  the  expansion  coefficients  C,  . .  satisfy  the  s tationary- state 

1 ,  n+fi  - 1 

Schrodinger  equation 


M 


)  ^^ii^ii'  “ii'  '  ^^ii' ^^i’ ,n+M-i'  ”  ° 


(5) 


i'-l 


The  matrix  elements  in  (5)  are  given  by 


i-1 

-  <i,n+M-i|H|  i,n+M-i>  -  (n+M-l))<n  +  ^ 

j-1 

^ii'  “ 


(5a) 
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A^,yn+m-i'  5^  +  A^,  _^yn+M- i ’+1  5^  ,  i  ^  i'  (5b) 


Here  we  have  defined  the  detuning  parameters 


=  hn  -  -  i  -  1.  2 . M-1 


(6) 


The  energy  eigenvector  1^  >  and  Che  corresponding  eigenvalue  E  are 

no  "  no 

determined  by 


H  <6  > 

'  no 


E 

no 


<6  > 
no 


(7a) 


> 


^i,n+M-l' 


i,n+M-i> 


(7b) 


The  subscript  o  in  (7a)  labels  the  eigenstates  in  the  subspace  in  question. 

If  the  atom  is  assumed  to  be  initially  in  the  M-th  level,  then  for  all 
time  t  >  0  the  system  can  only  be  found  in  subspaces  with  N  >  M.  There  are  M 
eigenstates  of  the  type  (7b)  in  every  subspace.  Hence  the  orthonormality 
condition  is 


I 


C 


o' 

i  ,n+M-i 


C 


o 

i,n+M-i 


o'o 


(8) 


and  the  completeness  relation  is 

®  M 

)  )  U  I  -  1  . 

[_  '  no  no ' 

n-0  0-1 


(9) 


We  can  now  solve  the  equation  of  motion  for  the  density  matrix. 


Its 


matrix  element  between  energy  eigenstates  satisfies  the  equation 


d  t  ,  n '  a  ' 


K  ^^nc7  ^n'ff' ^^na.n'a' 


(10) 


whose  solution  can  be  put  in  the  form 


•w  (E  -E  ,  ,)t 

na  n  a 


(11) 


The  initial  matrix  element  in  (11)  is  given  by 


(12) 


where  C  ^  is  the  inverse  matrix  of  the  transformation  defined  in  (4)  and  p 

nn’ 

is  the  initial  density  matrix  element  of  the  field. 

The  normalized  intensity  correlation  function  of  the  field  is  defined  as 


g^^^(t)  -  <a^^a^>/<a^a>^ 


(13) 


where 


<n>  -  <a^a>  -  Tr(pa^a) 


®  M 


M  M 


-S')  (n+M-i)  ^  M  ■  cf’  „  . 

l_  l_  i,n+M-i  i,n+M-i 


n-0  i-1 


a-1  o' -I 


X  p  , (0)  cos 

na ,  na 


(E  -E  ,)t 
na  na 


t2  2 
<a  a  > 


CO  M 

L  L 

n-O  i-1 


(n+M- i) (n+M- i- 1) 


I  (E-E)t 

\  \  /r,\  ng  na' 

)  )  C.  _  .  C.  ^  .  P  , (0)  cos  - 

/  i,n+M-i  i,n+M-i  na.na 

a-1  a'— 1 


.  (15) 


The  probability  of  finding  m  photons  in  the  field  at  time  t  is  given  by 


<i,n+M-i|p  i,n+M*i> 


n— m-M+l 


where  i  -  n+M-m  and  n  >  0.  Taking  the  inverse  transformation  of  (7b)  and 
plugging  it  into  (16),  we  find 


m  MM 

)  )  I 

n-m-M+1  j-1  j'-l 


-1  n+M-m  -1  n+M-m 


X  p  .  . , (0)  cos 

nj.nj' 


f(E  .-E  .Jt! 

_QJ _ nj___ 

K  J  ’ 


where  we  have  made  use  of  (11).  It  is  then  not  difficult  to  prove  that 


oo  M 


\  ^  I  'nj., 


(0)  -  1  . 


n-0  j-1 


The  statistical  properties  of  the  rield  are  completely  described  by  Eqs .  (13) 

and  (17),  which  we  shall  discuss  in  the  following  section. 


III.  Results  and  discussion 

Throughout  this  paper,  we  take  K  ”  1  in  our  numerical  calculation.  The 
above  formulas  are  derived  on  the  assumption  the  atom  is  initially  in  its 
highest  level  We  now  assume  further  that  the  field  is  initially  in  the 
coherent  state  with  mean  photon  number  n.  Thus  the  diagonal  matrix  element  of 
the  initial  density  matrix  of  the  field  is 


P 


nn 


— n 
n 

n ! 


-  n 
e 


(19) 


As  is  well  known,  when  Che  field  is  strictly  in  the  coherent  state,  the 

(2)  (2) 

equal- time  second-order  correlation  function  is  g  (t)  =  1.  When  g  (t)  < 

(2) 

1,  the  field  exhibits  anCibunching ,  and  when  g  (t)  >  1,  the  field  exhibits 

bunching.  The  time  evolution  of  the  correlation  function  calculated  from  (13) 

is  presented  in  Figs.  2  and  3  for  various  M  and  n.  Resonance  conditions 

A  -  0  have  been  assumed  here  for  simplicity.  It  is  observed  that  the 
i 

correlation  function  oscillates  irregularly  around  its  initial  value  of  1. 
Thus  the  field  exhibits  alternatively  bunching  and  antibunching  as  time 
develops.  For  a  fixed  M,  this  oscillation  amplitude  decreases  with  the 
increasing  n.  In  other  words,  the  stronger  the  initial  field  is,  the  weaker 
the  effect  of  bunching  or  antibunching  becomes.  Therefore,  strong  photon 
bunching  or  antibunching  results  if  the  initial  field  strength  weakens.  The 
same  conclusions  have  been  obtained  for  the  M  -  2  case . ^ 


10 


Physically,  it  is  not  difficult  to  understand  this  phenomenon.  Since  we 
have  assumed  that  the  field  is  in  the  coherent  state  initially,  at  any  later 
time  or  t  >  0,  the  total  field  is  the  sum  of  the  initial  field  plus  the 
radiation  field  of  the  stimulated  atom.  The  latter  acts  like  a  modulator  chat 
causes  the  initial  field  to  deviate  from  coherence.  The  smaller  the  n,  the 
easier  the  initial  field  is  modulated.  At  any  time  t  >  0,  the  total  field  is 
dominated  by  the  stimulated  radiation  field  which  is  an  oscillating  function 
of  time.  As  a  consequence,  the  total  field  oscillates  between  bunching  and 
antibunching .  If  we  compare  the  curves  for  the  same  n  but  different  M.  we 
find  that  the  effect  of  antibunching  decreases  when  the  number  of  atomic 
levels  increases.  Similar  results  have  also  been  found  for  a  common  upper  N- 

g 

level  atom  interacting  with  N-1  modes  of  cavity  fields. 

The  photon  number  fluctuation  is  directly  related  to  the  Q-parameter 
defined  as^^ 


<n2> 


-  <n> 


<n> 


-  1 


(20) 


When  Q  >  0,  the  fluctuation  is  large  and  the  photon  probability  distribution 
is  super  -  Poisson .  When  Q  <  0,  the  fluctuation  is  small  and  the  probability 
distribution  is  sub-Poisson.  The  lower  bound  of  Q  is  -1,  corresponding  to 
zero  fluctuation,  where  the  probability  has  its  sharpest  peak.  When  Q  -  0, 
the  photon  probability  has  a  Poisson  distribution,  indicating  a  coherent  state 
of  the  field.  The  time  evolution  of  the  Q-parameter  is  presented  in  Figs.  6 
and  7  for  different  M  and  n.  From  these  figures  we  can  see  that,  similar  to 
g(2)  ,  for  fixed  M  the  amplitude  of  the  oscillation  of  Q  becomes  smaller  when  n 
increases,  and  that  the  photon  probability  approaches  the  Poisson  distribution 


11 


when  n  increases.  For  fixed  n,  increasing  M  will  cause  the  photon  probability 

to  deviate  from  the  Poisson  distribution  and  tend  to  become  super- Poisson. 

Since  the  Q-parameter  is  related  to  the  normalized  Glauber  correlation 

(2) 

function  g  (0)  by 

Q  -  <n>[g^^^0)  -  1]  .  (21) 

it  is  clear  that  photons  exhibit  the  antibunching  phenomenon  whenever  the 

probability  distribution  is  sub-Poisson,  and  exhibits  bunching  phenomenon  when 

the  probability  has  a  super- Poisson  distribution.  It  should  be  noted  that  a 

(2) 

minimum  g  does  not  imply  a  minimum  Q  because  of  the  factor  <n>  in  (21)  , 

except  that  <n>  is  just  in  the  time  region  with  collapsed  oscillation.  This 

can  be  seen  by  comparing  Figs.  2,  4,  and  6  or  Figs.  3,  5  and  7  with  one 

another.  We  have  also  found  that  for  certain  n  and  M  the  field  can  be  kept  in 

(2) 

a  coherent  state  for  a  long  time  with  g  =1  and  Q  «  0  (see  Figs.  2(c)  and 
6(c)).  This  means  that  the  initially  coherent  field  can  recover  its  coherence 
after  the  interaction  with  the  atom  takes  place.  This  is  in  agreement  with 
the  conclusions  of  Ref.  6. 

Next  we  look  at  the  case  in  which  thg  detunings  are  nonzero.  We  assume 
that  |a^|  ^  0  (i  -  1,  2,  ...,  M-1)  but  ^  -  0 .  That  is,  multiphoton 

transitions  involving  less  than  M-1  photons  are  all  detuned,  and  only  the 
(M-l)-fold  dfc  -.lerate  photon  transition  satisfies  the  resonance  condition. 
The  correl  - ^  ...v^n  function  is  calculated  for  M  -  3  and  9,  n  -  3  and  [a^I  -  0.5 
and  1.5.  The  ..--vults  are  plotted  in  Figs.  8  and  9.  We  see  from  these  figures 
that  for  a  given  M,  a  larger  detuning  results  in  weaker  antibunching  but 
stronger  bunching.  This  feature  becomes  even  stronger  as  M  becomes  larger. 
On  the  other  hand,  when  the  detunings  Ia^^I  are  fixed,  then  larger  M 
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corresponds  to  weaker  antibunching.  This  means  that  the  more  photons  that  are 
involved  in  the  degenerate  multiphoton  transition,  the  weaker  the  antibunching 
becomes  in  the  process. 

In  addition,  we  have  also  examined  the  sub-Poisson  distributions  of  the 
photons  with  the  same  M  but  different  n,  where  we  find  sharper  peaks  in  the 
stribution  when  the  inititai  field  strength  is  weaker  or  when  n  is  smaller, 
-'ted,  however,  that  n  should  be  kept  above  a  certain  value,  say,  >  3, 
-therwise  the  distribution  approaches  the  shape  of  a  thermal 
distri..  i.  We  have  also  computed  the  photon  distribution  for  nonresonant 
cases  (numerical  results  are  not  present  d  here),  where  we  have  found  that  the 
detunings  tend  to  weaken  the  sub-Poisson  and  super- Poisson  effects. 

IV.  Conclusions 

We  have  shown  that  multiphoton  processes  lead  to  a  minimum  fluctuation 

in  the  field  intensity.  Although  such  a  markedly  reduced  fluctuation  occurs 

only  at  some  particular  time,  it  can  be  achieved  by  adjusting  the  velocity  of 

the  atom  injected  into  the  optical  cavity,  because  the  JC  model  is  realized  in 

9 

practice  by  the  micromaser.  In  other  words,  the  most  stable  intensity  of  the 

output  field  can  be  obtained  by  adjusting  the  duration  that  the  atom  stays 

12 

within  the  cavity.  This  is  also  true  for  the  squeezing  of  the  field. 

Acknowledgments 

This  research  was  partially  supported  by  the  Office  of  Naval  Research, 
the  National  Science  Foundation  under  Grant  CHE-8620274,  and  the  Air  Force 
Office  of  Scientific  Research  (AFSC) ,  United  States  Air  Force,  under  Contract 
F49620-86-C-0009 .  The  United  States  Government  is  authorized  to  copy  and 


13 


distribute  reprints  for  governmental  purposes  notwithstanding  any  copyright 
notation  hereon. 


1 


14 


References 


1. 

E. 

T. 

Jaynes 

and  F.  W.  Cummings,  Proc.  IEEE  51 . 

89  (1963). 

2. 

F. 

W. 

CununiT 

hys.  Rev.  A  140,  1051  (1965). 

3  . 

P. 

Meystre , 

E.  ueneux,  A.  Quattropaniand  and  A. 

Faist,  Nuovo  Cimento  25B . 

521  (1975) . 

4.  R.  R.  Puri  and  G.  S.  Agarwal,  Phys.  Rev.  A  23,  3610  (1986). 

5.  M.  Hillery,  Phys.  Rev.  A  35,.  4186  (1987). 

6.  Z.  D.  Liu,  X.  S.  Li  and  D.  L.  Lin,  Phys.  Rev.  A  36,  5220  (1987). 

7.  D.  L.  Lin,  X.  S.  Li  and  Y.  N.  Peng,  Phys.  Rev.  A  39,  1933  (1989). 

8.  M.  Kozierowski  and  A.  S.  Shumovsky,  Physica  145A.  290  (1987). 

9.  G.  Rempe,  H.  Walther  and  N.  Klein,  Phys.  Rev.  Lett.  52.,  353  (1987);  D. 

Meshhede,  H.  Walther  and  G.  Mailer,  Phys.  Rev.  Lett.  54.  551  (1985). 

10.  M.  Brune,  J.  M.  Raimond,  P.  Goy,  L.  Davidovich  and  S.  Haroche ,  Phys.  Rev. 

Lett.  59,  1899  (1987). 

11.  F.  L.  Li,  D.  L.  Lin,  T.  F.  George  and  X.  S.  Li,  Phys.  Rev.  A,  in  press. 
This  paper  is  hereafter  referred  to  as  I. 

12.  F.  L.  Li,  X.  S.  Li,  D.  L.  Lin  and  T.  F.  George,  J.  Phys.  B,  in  press. 

13.  L.  Mandel,  Opt.  Lett.  4,  205  (1979). 


15 


Figure  captions 

1.  Schematic  diagram  of  the  atomic  levels. 

(2) 

2.  Time  evolution  of  the  correlation  function  g  under  resonance  conditions 

=  0  for  M  -  3:  (a)  n  -  3,  (b)  n  -  5,  (c)  n  -  10. 

3.  Same  as  Fig.  2  except  M  =  9. 

4.  Time  evolution  of  the  mean  photon  number  <n>  for  -  0  and  M  -  3:  (a)  n 

-  3,  (b)  n  -  5,  (c)  n  -  10. 

5.  Same  as  Fig.  4,  except  M  -  9. 

6.  Time  evolution  of  Handel's  Q-parameter  for  A^  -  0  and  M  -  3:  (a)  n  -  3, 

(b)  n  =  5,  (c)  n  =■  10. 

7.  Same  as  Fig.  6,  except  M  -  9. 

(2)  - 

8.  Time  evolution  of  g  for  M  -  3  and  n  -  3.  The  resonance  conditions  are 

not  satisfied  by  the  individual  transitions,  but  -  0,  where 

.  A^  >  0,  .  <  0  (i  -  1 . M-1):  (a)  |a^|  -  0.5,  (b) 

2  ~ 

|aJ  -  1.5. 

9.  Same  as  Fig.  8  except  M  =  9. 
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